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ABSTRACT 

One goal of CMB data analysis is to combine data at different frequencies, angular 
resolutions, and noise levels in order to best extract the component with a Plankian 
spectral behaviour. A multi-frequency Wiener filtering method has been proposed in 
this context by Bouchet, Gispert and Puget (1995) and in parallel by Tegmark and 
Efstathiou (1996). As shown in Bouchet and Gispert (1998a), this linear method is also 
convenient to estimate a priori, given a sky model and an experimental description, 
the residual errors on the CMB power spectrum assuming the foregrounds have been 
removed with this method. In this paper, we extend the method to the case when 
additional polarisation data is available. In particular, we derive the errors on the 
power spectra involving polarisation and show numerical results for the specifications 
of the future CMB space missions MAP and Planck * when it is assumed that the 
Galactic synchrotron and dust emission are respectively about 40% and 10 % polarised. 
We consider two underlying models for our study: we take a standard CDM model 
with t = 0.1 for the extraction of E-mode polarisation and ET cross-correlation ; for 
-B-mode polarisation we consider a tilted CDM model with n s = 0.9, tit = —0.1 and 
T/S = 0.7. We find that: (1) The resulting fractional errors on E mode polarisation 
and TE cross-correlation power spectra are ^ 10-30 % for 50 t & 1000 for Planck. 
The fractional errors are between 50% to 150% for I < 50, (2) The corresponding 
fractional errors for MAP are > 300 % for most of the t range, (3) the Wiener filtering 
give extraction errors < 2 times the expected performance for the combined sensitivity 
of all the channels of Planck . For MAP, the corresponding degradation is ~ 4. (4) 
if, instead of individual modes, one considers band-power estimates with a logarithmic 
interval Al/£ = 0.2 then the fractional error for MAP drops to 100% at the Doppler 
peak around I ~ 300 for the ET signal, and (5) The fractional error for B-mode 
polarisation detection is ^ 100% with Planck for I < 100. A band-power estimate 
with Aijl = 0.2 reduces the fractional errors to <> 25% for 20 < I < 100. 



1 INTRODUCTION 



Ever since the detection of CMB temperature anisotropies 
by the DMR experiment aboard the COBE satellite at an- 
gular scales > 7° (Smoot et al. 1992, Bennett et al. 1996), 
there has been a surge of activity on both theoretical and 
experimental front. Several detections at smaller angular 
scales have been reported (for a recent compendium, see 
Lineweaver and Barbosa 1998). There is great interest in 
the upcoming satellite projects MAP and Planck. These 
all-sky experiments, it is hoped, will determine the angu- 
lar power spectrum of CMB temperature fluctuations at all 
scales greater than a few arc minutes. Various theoretical 
analyses have shown under fairly general assumptions that 
this could determine cosmological parameters, Q, Qb, h, etc, 
with unprecedented precision (Jungman et al. 1996, Zal- 
darriaga et al. 1997, Bond et al. 1997). In addition, there is 
growing optimism that these satellite projects might detect 



the very small expected signal from the polarised component 
of the CMB. The information from CMB polarisation could 
then be combined with temperature signal 1) to check the 
self-consistency of the predictions of the underlying theoreti- 
cal model, 2) to further break the degeneracy between some 
cosmological parameters and 3) to help in unambiguously 
detecting a tensor-induced component of CMB fluctuations 
(Seljak & Zaldarriaga 1997, Seljak 1997, Kamionkowski & 
Kosowski 1997). 

A major hurdle in extracting the primary CMB signal 
from data, apart from noise, is the presence of poorly known 
Galactic and extra-galactic foregrounds. However, as the 
foregrounds differ from the CMB emission in both frequency 
dependence and spatial distribution, one can reduce their 
level by proper combination of the multi-frequency data of 
a CMB experiment. Bouchet et al. (1995) and Tegmark & 
Efstathiou (1996) proposed a particular linear scheme, based 
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on the traditional Wiener filtering method to take full ad- 
vantage of this fact. The performance of the method has 
been assessed through detailed numerical simulations per- 
formed in the context of the Planck preparation (Gispert 
& Bouchet 1996, Bouchet & Gispert 1998b). It was shown 
that the residual contamination after cleaning the map is 
much smaller than the CMB primary signal, and therefore 
the foregrounds may not be a major obstacle in the extrac- 
tion of CMB temperature angular power spectrum. 

The CMB polarisation signal is expected to be one to 
two order of magnitudes below the temperature signal. And 
it is likely to be comparable to the achievable experimen- 
tal noise in the current experiments. The presence of fore- 
grounds should, of course, make this detection even more 
difficult. In our Galaxy, the synchrotron emission is highly 
polarised and would constitute a major foreground for mea- 
surements made in the Rayleigh- Jeans part of the CMB 
spectrum (e.g. MAP and Low-Frequency Instrument -LFI- 
of Planck). The Galactic dust emission is also observed to 
be polarised and it is the dominant foreground in the Wien 
part of the CMB spectrum (case of High- Frequency Instru- 
ment -HFI- of Planck). Prunet et al. (1998) modelled and 
estimated the level of dust polarised emission at high Galac- 
tic latitudes, and compared it to the polarised component 
of the CMB signal. They showed that the scalar-induced 
polarisation (E-mode) power spectrum and temperature- 
polarisation cross-correlation are likely to be much larger 
than the dust polarised emission at high Galactic latitudes. 
However, the tensor-induced (B-mode) signal is at best com- 
parable to the foreground contamination level. The extrac- 
tion of such a signal in the presence of comparable instru- 
mental noise, even with small foregrounds, is trickier than 
the corresponding temperature case, where the signal is 
much bigger than the instrumental noise. 

In this paper, we extend the mult i- frequency Wiener fil- 
tering method to include the polarisation and temperature- 
polarisation cross-correlation. The goal of this exercise is to 
quantify errors in estimating various power spectra when 
Galactic polarised foregrounds are present. We describe the 
method in detail in next section. In § 3, we give a description 
of the polarised foregrounds emission which is used in § 4 
to give the resultant errors on the power spectra estimates 
in the case of MAP and Planck. We summarize our results 
and discuss their limitations and applicability in § 5. 



2 MULTI-FREQUENCY WIENER FILTERING 
ON CMB DATA 

We present here an extension of the multi-frequency Wiener 
filtering technique applied earlier on CMB temperature data 
(Bouchet et al. 1995, Tegmark & Efstathiou 1996, Bouchet 
& Gispert 1998a, hereafter BG98). The goal of this extension 
is to see the effect of including polarisation data. 

Let us denote the observed data at different frequencies 
as, yl, where v indicates the frequency of the instrumental 
channel, and i the nature of the observed field (tempera- 
ture or .E-mode polarisation. The B mode polarisation has 
a vanishing cross-correlation with both T and E and there- 
fore can be treated separately). The E and B modes are 
the 'divergence' and 'curl' part of the polarisation tensor 
and are linear combinations of the usual Stoke's parameters 



Q and U. The advantage of using these variables is that 
the B mode polarisation vanishes for scalar-induced pertur- 
bations (for details see Seljak 1997). The data points yl, 
for a given i, take contributions from various Galactic and 
extra-galactic sources, apart from the primary CMB signal 
and instrumental noise. Throughout this paper, we adopt 
the convention of Tegmark and Efstathiou (1996) and take 
yl to mean the 'observed' temperature fluctuation, i.e., the 
monopole part is first subtracted from the observed surface 
brightness, which is then divided by dBo/dTo, Bo being the 
surface brightness of CMB at T — Tq. Let us call x 3 p the 
contribution of the field j due to process p, this is the quan- 
tity we want to recover from the observational data yl. We 
assume there is a linear relation between them: 



i a ij j , ii 

— A vp * X p + V , 



(1) 



where A]? p is the instrument response kernel, and bl is the 
detector noise level per pixel for the full mission time (all 
repeated indices are summed over). As discussed in BG98, 
it is easier to deal with the spherical harmonics transform of 
Eq. (Q), as this non-local equation in pixel space transforms 
into a algebraic linear equation in multipole space. From 
now on, all spatially dependent variables are then expressed 
in multipole space. Eq. {jlh translates into: 



yl(l,m) = A l J p (l,m)x P (l,m) + bl(l,m) 



(2) 



The problem now is to construct an optimal estimator xp. 
We chose this estimator to be linear in the observed data: 

K = W;iyl. (3) 

The reconstruction error for a given field and process (ep) 2 
:./'.. — x p ) 2 \ can be written as: 

( £ J>) 2 = { W P» A l k p> ~ $PP' $ik) {Wpu'A^p" ~ S pp "S im ) (x p , 

+w;:w;t,(b a X,); 



Wpi is chosen so as to make the the variance of the recon- 
struction error minimal. The derivatives of the error with 
respect to the filters coefficients Wp^ should then be zero. 
This condition can be expressed as: 

[(SjcSpuAl*,) (Wp V ,A l ™p„ ~8pp„5 im ) 

+ (S^'SicA 1 ™^,) (WplAl h p , -SppiSik)] {xp'Zp} 

+ (S^SaaWfr + S^ScbWg) {b a X'} = 0. (5) 
Rearranging terms we get the final equation on the filters: 
^''Wp^l7 p »(VV) + Wp„,(bp,b b „,} = Afp,(x k v ,x p ). (6) 



2.1 Implementation of foregrounds removal 

Eq. (|^) is valid for the general case in which various pro- 
cesses, fields, and corresponding instrumental noises could 
be correlated. We consider here only uncorrelated processes 
and uncorrelated noises between different fields and chan- 
nels. We allow for the correlation between the two fields T 
and E. These conditions can be expressed as: 



/ l l i 

\Xp'Xp!l j 

I 1 2 \ 
\Xp/Xp// j 

I 2 2 v 
\XpiXpii j 



Op' p" Op/ iAp/p// 
c fiTE _ A.TE 

Op' p" Op/ iip/p// 

e s~lE * E 

Op' p" Op/ iAp/p// 



(7) 
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<«V&*> = (8) 

With these conditions, Eq. (^) can be written as a system 
of four matrix equations: 

W 11 (A n A T tr(A n ) + N T ) 



+ W 12 A 22 A Jfi tr(A 



W 12 (A 22 A E tr(A 2 



+ 



W 11 A 11 A TB tr(A 2 



1 A ii A J 

r21 Mun 



: A T tr(A n ) 
-N B ) 

: A TB tr(A 22 



(A A tr(A ) + N j 
+ W 22 A 22 A TB tr(A 11 ) = A TB tr(A 11 ) 

W 22 (A 22 A E tr(A 22 )+N B ) 
+ W 21 A 11 A TB tr(A 22 ) = A B tr(A 22 ) 



(9) 



(10) 



(11) 



(12) 



These equations can be solved by substitution. One can 
readily verify that the equation for W 11 reduces to the form 
earlier derived by Bouchet et al. (1995) and Tegmark & Efs- 
tathiou (1996) when the cross-correlation between the fields 
is switched off. 

Bouchet et al. (1996) defined a quantity they termed 
'quality factor' to quantify the merit of extraction of the 
signal corresponding to a given process. A straightforward 
generalization of this quality factor, valid for multiple fields, 
can be written as: 

QZ, = p^ = w;U k J P >'(x l P ^ p ,) (i3) 

(x P x p >) 

where we have used Eq. (|^) to write the second equality. 
Since this is is just the ratio of the spectra of the minimum- 
variance estimated map to the spectra of the real one, it can 
be viewed as the effective window function of the experi- 
ment. 

Eq. jl3| can be expanded to yield: 



(A T )- X (W 



11 A 11 A T + W 12 A 22 A TB ) 



(A E ) \ W 21 A' A ; 



,■21 i 11 a TE 



W 22 A 22 A B ) 



(A TB ) _1 ( W n A n A TB + W 12 A 22 A B ) 



(14) 
(15) 
(16) 
(17) 



Q 

Q 22 
Q 12 

Q 21 = Q 1 

Q 11 and Q 22 can readily be interpreted as the quality of 
the reconstruction of temperature and polarisation maps, re- 
spectively. As expected in the presence of cross-correlations, 
the quality factor of either field is better than in the case 
without cross- correlations. Although the reconstruction of 
temperature maps is only slightly changed by the cross- 
correlation term (the term proportional to W 12 in Eq. (|il|)), 
the quality of the polarisation reconstruction gets a big 
boost from the presence of temperature-polarisation cross- 
correlation and almost half the contribution to Q 22 comes 
from the W 21 term in Eq. @. However the meaning of the 
term Q 21 (and Q 12 ) is not apparent. Much of the contribu- 
tion to Q 12 comes from the term with W 11 , and therefore 
it is very close to the quality factor for the extraction of 
temperature and is nearly independent of the polarisation 
noise. It is not surprising as it merely tells us that to opti- 
mally reconstruct the cross-correlation one needs to throw 
out the noisy data, i.e. the polarisation information! 



However, the quantity of interest is the error in the 
extraction of the power spectrum of cross-correlation which 
should not be confused with Q 12 . To get a real idea of the 
error bars of the different spectra, we must define estimators 
of those power spectra from the filtered data, and compute 
their covariances. As we shall see, while Q 11 and Q 22 directly 
give the covariance of the E and T power spectra, a more 
complicated expression is needed for the covariance of the 
ET power spectrum. 

2.2 Unbiased estimators of power spectra 

Eq. (|^) is the data obtained after performing Wiener filtering 
on the multi-frequency maps. Our aim in this section is to 
use this data to write an unbiased estimator of the power 
spectra, which is defined such that {C p ) = (x p x p ). From 
Eqs. (Q) and (^), the power spectrum of x p can be written 
as: 



(44) = w;lwp [Al n p ,A™i„(x;xD + (&i^>] , 

which can be expressed as: 

\XpXpi \~p * Op ) ^°lp ^p 

K^p ^p i ^p ) ^p 



(18) 



(44> = (z?cj E + b?) = , 



(19) 



The meaning of the first equality in Eqs. ([[9]) can be easily 
understood: the Wiener-filtered power spectrum for a given 
process can be expressed in terms of the true underlying 
power spectrum smeared by foregrounds {Z % p terms) plus 
the noise. However, unlike the case with no foregrounds, 
The 'noise' terms b z p have contribution not only from the 
instrumental noise but also from residual foregrounds from 
other processes and fields. The second equality comes from 
Eq. (fL3|) . These equations can then be used to write unbiased 
estimators of the power spectra: 



d * = z]? ( W+T ^ H*pM4( r 

\ m 

d ? = ~z\* ( 2T+1 Pp( m ) £ p( r 

\ m 



1 



\Xp(rn)x p (m)\\ — b. 



(20) 
(21) 
(22) 



From the unbiased estimators, one can readily compute the 
covariances of the various power spectra: 



Cov((7 p T ) = ^—^{clQflZf) 2 (23) 
^Tj-j- (C p Qf /Z p 2 ) 2 (24) 

1 (iQl'ncD 2 + Ql 1 Ql 2 CpC E ) 



Cov(C p 



Cov(c; ' 



21- 



4*5) 



Given the instrumental noise and the expected level of fore- 
grounds, Eqs. p3[), (p^), and ( pB[ ) can be used to estimate 
the precision with which the power spectra of various quan- 
tities can be determined. In the next section, we use the 
specifications of future experiments MAP and Planck to 
estimate the fractional errors on the power spectra. In calcu- 
lating the covariances, we have assumed all the processes — 
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CMB and foregrounds — to be Gaussian. If a fraction of sky, 
/sky, is covered then the corresponding expressions for the 
covariances can be obtained by dividing the equations (p3|), 
(pi|), and p5[ ) by / s k y - Throughout this paper we assume 
complete sky coverage, ie. / s k y = 1. 



3 POLARISED FOREGROUNDS 

The most dominant Galactic polarised foregrounds are ex- 
pected to be the polarised components of synchrotron and 
dust. In addition, it is possible that the free- free emission 
is also polarised at 10% level, which could be a further de- 
terrent to extracting the CMB polarisation. We neglect the 
possibility of polarised free-free emission in this paper (for 
more details see Keating et al. 1997). 

The galactic synchrotron radiation originates from the 
interaction of cosmic ray particles with the galactic magnetic 
field and is known to dominate the galactic radio emission 
for frequencies < 10 GHz. In theory, this radiation can be 
— 70 % polarised, for the observed energy spectrum of the 
cosmic ray particles (Rybicki & Lightman 1979). However, 
because of the cancellation due to incoherent addition of po- 
larised component along any line of sight, Faraday depolar- 
isation, and non-uniform magnetic fields, the true polarised 
percentage is observed to be ~ 40% with a dependence on 
galactic latitude (Burn 1966; Spoestra 1984). The galactic 
emission in FIR and millimeter wavelengths is dominated 
by dust emission. The dust particles align themselves with 
the interstellar magnetic, and because the dust particles are 
not spherical the resulting FIR and millimeter emission is 
polarised (Hildebrand & Dragovan 1995). 

Polarised dust : The Galactic dust is seen to be polarised 
at a small level (^ 10%). This is likely to be the major fore- 
ground at frequencies at which the Planck HFI will op- 
erate. Though there exist no observations of polarised dust 
at high Galactic latitudes, it is possible to model this emis- 
sion using observations at smaller Galactic latitudes. Prunet 
et al. (1998) modelled this distribution and computed the 
power spectra of polarisation and temperature-polarisation 
cross-correlation. They showed that the relevant power spec- 
tra can be fitted as: 



C d E St {£) 


= 8.9xl0~ 4 r 


13 ( M K) 2 


(26) 


C B ust (£) 


= i.oxio~ 3 r 


- 1 ' 4 (mK) 2 


(27) 


Cet{1) 


= i.7xio~ 2 r 


■ 1M (/,K) 2 . 


(28) 



These power spectra are normalized at 100 GHz. They cor- 
respond to Galactic latitudes between 30° and 45° and are 
taken here as representative of the all sky average (for more 
details on this and other related issues, see Prunet et al. 
1998). The dust emissivity is assumed to be proportional to 
v 2 with a temperature of 17.5 K (Boulanger et al. 1996). For 
comparison, we recall that the temperature power spectrum 
at 100 GHz can be fitted by (BG98): 



176 



( M K) 2 



(29) 



Polarised Synchrotron emission: This foreground will 
greatly undermine the performance of MAP and Planck 
LFI in the detection of CMB polarisation. The existing 
maps of polarised component of Galactic synchrotron shows 
this emission to be polarised at a level between 20 to 60%, 



depending on the Galactic latitude, at radio frequencies 
< 1.4 GHz (Spoestra 1984). These can be extrapolated to 
millimeter wavelengths using the frequency dependence of 
synchrotron emission (Lubin & Smoot 1981). However, such 
a procedure is fraught with uncertainties: (1) The syn- 
chrotron spectrum is not well known up to millimeter fre- 
quencies (for recent attempts to measure the spectrum up 
to 10 GHz see Platania et al. 1997), (2) At low frequen- 
cies (< 1.4 GHz) the polarised emission is not optically thin 
because of Faraday depolarisation. This means that the ob- 
served emission suffers from substantial depolarisation (and 
also the spatial distribution is affected because of the rota- 
tion of the polarisation axis). The Faraday depolarisation 
is proportional to v~ 2 and becomes negligible at millimeter 
wavelengths. Therefore, the degree of polarisation at these 
frequencies is expected to be higher; in addition it will have 
a different spatial distribution. So the existing data needs 
to be corrected for these effects when an extrapolation to 
higher frequency is performed. It is not easy to do so with 
the present data. For the lack of data, we assume the syn- 
chrotron emission to be polarised at 44% level with the same 
spatial distribution (^-dependence) as the unpolarised emis- 
sion. We assume both E and B mode power spectra to cor- 
respond to this level of polarisation. It is probably justified 
because the modelling of dust polarised emission also shows 
comparable emission for these mode. Furthermore, we as- 
sume perfect cross-correlation between _E-mode polarisation 
and temperature. This is likely to be the case because both 
the polarised and unpolarised emission depend on the square 
of the magnetic field (its component in the plane of the sky) 
and the emission at high Galactic latitudes is mostly domi- 
nated by one structure. We obtain: 



cum 

where 



o.2 x c s T yn (£) 

0.2 x C s T yn (£) 
0.44 x C 3 ^ n {£) 



C s T yn (£) =4.5r J (nK) 2 



(30) 
(31) 
(32) 

(33) 



at 100 GHz (BG98). One of the aims of taking high 
level of synchrotron polarised emission and perfect cross- 
correlations is to consider the 'worst possible' case for the 
performance of MAP and LFI. This case should be con- 
trasted with the CMB case in which the cross-correlation 
is ~ 1/3 of the perfect cross-correlation. And therefore the 
assumed perfect cross-correlation in any foreground means 
that the total signal is biased in favour of the foreground and 
would make the extraction of CMB cross-correlation more 
difficult. However, an important goal is also to ask how well 
can these foregrounds be extracted using the same experi- 
ment that attempts to measure the CMB polarisation. We 
address this question in a later section. As mentioned above, 
the free-free emission can also be polarised at a small level 
but for the assumed level of synchrotron emission it will be 
sub-dominant to it at all frequencies. 

In Fig. |lj we compare the expected CMB E and ET sig- 
nal with the assumed level of foregrounds. The CMB signal 
dominates the foregrounds for frequency channels between 
44 GHz and 217 GHz. The lower frequency channels of MAP 
(22 and 30 GHz) and LFI (30 GHz) and the higher frequency 
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channels of HFI (545 GHz) will be dominated by polarised 
synchrotron and dust, respectively; these channels will help 
in an accurate determination of these foregrounds. 

Recently, it was argued by Draine and Lazarian (1998) 
that a part of galactic foreground between 10 and 100 GHz 
could be contributed by non-thermal emission from rotating 
dust grains. If this emission is polarized at the same level as 
we assume in this paper (Eq. p8[ ) then it could add to the 
polarised galactic emission. However, this emission is sub- 
dominant to the polarised synchrotron emission we assume 
in this paper. And therefore we neglect the effect of spinning 
dust particles in our discussion in this paper. 

Apart from the polarised foreground and in addition to 
the Galactic unpolarised foregrounds — dust, free-free, and 
synchrotron — we include several other extragalactic un- 
polarised foregrounds in our study: The thermal Sunayev- 
Zeldovitch effect from clusters of galaxies, the infra-red point 
sources, and the radio point sources (for details see Bouchet 
& Gispert 1998a). We do not consider any extragalactic po- 
larised foregrounds. 



4 EXTRACTION OF POLARISED SIGNAL 

Quality factors: As mentioned in the last section, a relevant 
quantity for determining the merit of signal extraction is 
the 'quality factor' (Eq. (|l3[)). The quality factors for tem- 
perature and i?-mode polarisation are shown in Figs. (^|) 
and (^) for various processes. It should be noted that in 
these figures we only show the quality factor for unpolarised 
synchrotron and dust though we include several other unpo- 
larised Galactic and extragalactic foregrounds in our anal- 
ysis (see §3). Throughout this and the next subsection, we 
take the underlying models for calculating various power 
spectra to be a variant of sCDM model with n s — 1, the ra- 
tio of scalar to tensor quadruple T/S = (no tensor signal), 
and the optical depth the last scattering surface r = 0.1. 
Finite optical depth to the last scattering surface enhances 
the polarisation signal for I ^ 10, a part of the spectrum 
which is of great interest for breaking degeneracies between 
various cosmological parameters (for details see Zaldarriaga, 
Spergel, & Seljak 1997). Throughout this paper, we use the 
CMB Boltzmann code CMBFAST to generate CMB fluctu- 
ations (Seljak & Zaldarriaga 1996). 

As expected, the temperature signal can be extracted 
far more cleanly than the polarised signal. The quality fac- 
tor for the extraction of CMB drops exponentially as the 
£ approaches the effective beam width of the respective ex- 
periment. Another noticeable feature in temperature quality 
factors for various experiments is that the spatial distribu- 
tion of Galactic dust emission can be discerned almost as 
well as the CMB signal using Planck HFI. This is largely 
owing to the presence of polarised channel at 545 GHz. The 
signal at these channels will be dominated by Galactic dust 
emission; and they have sufficiently low noise levels and high 
enough angular resolution to allow a good determination of 
the power spectra of the dust emission up to £ ~ 1000. This 
should be contrasted with the extraction of unpolarised syn- 
chrotron emission, which is one of the dominant contaminant 
for v < 90 GHz. Unfortunately neither MAP nor LFI can ex- 
tract this signal well. It is because (a) it does not dominate 
the signal at any frequency of either MAP or LFI (b) the 



signal is the strongest is the lowest resolution channels, and 
(c) the unpolarised synchrotron and Galactic free-free emis- 
sion have very similar spatial distribution (assumed to be 
the same in this paper) and nearly the same frequency de- 
pendence, which makes it difficult to extract either of them. 

The oscillating nature of quality factors for CMB polar- 
isation attests to the fact that the signal is barely above the 
noise level. Another important thing to notice in the figures 
is the quality of extraction of processes like the polarised 
component of dust and synchrotron. For HFI, the polarised 
dust can be extracted better than the CMB for much of 
the ^-range. As in the temperature case, it is possible be- 
cause of the presence of a polarised channel at 545 GHz. 
Equivalently, the lowest frequency channels of MAP and 
LFI serve as templates for polarised synchrotron. It should 
be noted that the synchrotron polarised component can be 
extracted much better than its unpolarised counterpart, be- 
cause, whereas the former dominates all other signals at the 
lowest frequency channels of MAP and LFI, the latter, as 
discussed above, doesn't dominate the total signal at any 
frequency, at least in our modeling. 

We shall see below how the information on quality fac- 
tor translates into errors on power spectra for various pro- 
cesses. 

4.1 Errors on power spectra 

We use Eqs. (p3|), (p4|), and (E9) to estimate the fractional 
errors in the extraction of the CMB _B-mode and TE cross- 
correlation power spectra. The results are shown in Figs.(Q), 
(g), and ^ for the specification of various experiments. For 
comparison, we also plot the expected errors using the best 
channel (90 GHz channel for MAP, 100 GHz channel for LFI, 
and 143 GHz channel for HFI and full Planck) and the 
combined sensitivity of all the channels of each experiment 
(Bond et al. 1997). 

For temperature signal, all the experiments give similar 
results for £ < 300. It is expected as the signal of temper- 
ature fluctuation is so much above the noise level for all 
the experiments in this ^-range that additional sensitivity 
does not lead to additional precision in power spectrum es- 
timation. The only source of error in this ^-range is cosmic 
variance which is obviously independent of the experiment. 
This information could also be gleaned from the quality fac- 
tor plots. For £ > 300, the extraction of temperature signal 
depends on the relative beam widths of relevant channels 
of various experiments. As expected, HFI performs best be- 
cause it will have channels with 5' resolution; it is followed 
in performance, in that order, by LFI and MAP. It should 
be noted that the results using Wiener filtering generally 
lie between the performances of best channel and combined 
sensitivity. 

Unlike the temperature fluctuations, the extraction of 
polarisation and temperature-polarisation cross-correlation 
depends sensitively on the pixel noise, as is evident from 
comparison between fractional errors on these quantities for 
various experiments. (The sharp spikes in ET plots is not 
a property of the extraction error but merely indicate that 
the signal vanishes at these values of I in the underlying 
theoretical model.) The performances of LFI and HFI are 
similar for £ 600 largely because of comparable polarised 
sensitivity of the 100 GHz LFI and 143 GHz HFI channels. 
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For larger £, HFI performs better because of its higher angu- 
lar resolution. Both these experiments should extract both 
the E-mode power spectrum and ET cross- correlat io n to 
20-30% precision, i.e. signal-to-noise, [Cnj W Cov(Cf)), of 
4-5, for 50 ^ £ 1000. The signal-to-noise is much smaller 
at larger scales, apart from cosmic variance, because the 
CMB signal becomes comparable to both foregrounds and 
pixel noise. A comparison with the results from the com- 
bined sensitivity of all channel shows that the presence of 
foregrounds degrade the signal extraction by a factor < 2 
for much of the £ range. 

From the figures, it is clear that MAP is hampered 
not only by foregrounds but also by its sensitivity. If the 
foregrounds were neglected the combined sensitivity of all 
its channels might enable a marginal detection of this sig- 
nal. However, our analysis shows that the presence of fore- 
grounds will degrade this detection by a factor ~ 4, and the 
resultant fractional errors on E and ET cross-correlation 
power spectra will be > 300% in most of Grange except 
I < 10. And therefore it seems unlikely that MAP could 
give a positive detection of either E and ET signal for indi- 
vidual modes. 

Band-Power estimates : A way to reduce the 'noise' fur- 
ther is by considering band-power estimates around a given 
£ (see Bond 1996 for relevant definitions). Band power taken 
over a logarithmic interval of A£/£ results in the reduction of 
errors by a factor of \J~H£ around any £. In Figs (Q) and (||), 
expected 1<t measurements are plotted for various experi- 
ments, with band powers taken over a logarithmic interval 
of 20%, i.e. A£/£ = 0.2 for Planck HFI and LFI. As is 
clearly seen in the pictorial representation of these figures, 
LFI and HFI should make a fairly accurate detection of E 
and TE power spectra for £ ^ 750. For MAP these estimates 
show that the errors for ET signal drop to ^ 100% near the 
Doppler peak at i ~ 300. 

It should be borne in mind that in the absence of any 
data on the power spectrum of polarised synchrotron we 
took the synchrotron power spectrum to have the same 
shape as the temperature power spectrum with 44% polari- 
sation and a perfect cross correlation between E and T . By 
reducing the level of synchrotron foreground and the ET 
cross-correlation it should be possible to get lower errors. 
On the other hand, a higher level of polarised synchrotron 
should allow a better extraction of synchrotron itself, which 
might be subtracted using methods other than Wiener filter- 
ing. Therefore, there is trade-off between a small and large 
assumed level of foregrounds: a sub-dominant foreground 
will give smaller formal errors on the extraction of CMB po- 
larisation but these foregrounds themselves will be elusive, 
thereby making it harder to quantify the errors; on the other 
hand, a larger foreground can be detected and probably sub- 
tracted more efficiently using, for instance, its non-Gaussian 
nature. These considerations makes it worthwhile to quan- 
tify the errors on foreground extraction. 

4.2 Extraction of polarised foregrounds 

One of the important goals of future mult i- frequency experi- 
ments is to detect and subtract the foregrounds which hinder 
the determination of CMB signal. At present one has to ex- 
trapolate spectral information on various foregrounds from 
radio and FIR wavelengths to millimeter and sub-millimeter 



wavelengths at which the CMB experiments operate. A sim- 
ilar extrapolation is required on the spatial distribution of 
foregrounds. For instance, the synchrotron spectrum is de- 
termined only for v ^ 10 GHz (Platania et al. 1997) while 
the spatial information is known for angular scales > 0.5° 
from the existing all-sky maps (Haslam et al. 1981). While 
dust maps are available with an angular resolution of 5' one 
needs to extrapolate the dust spectrum from 60 fim to mil- 
limeter wavelengths for a comparison with the CMB signal 
(Neugebauer et al. 1984). It has been pointed out that such 
an extrapolation, though useful, might result in large errors 
in a high-sensitivity CMB experiment (Brandt et al. 1994). 
Therefore, it is of paramount importance to determine the 
spectrum and spatial distribution of the foregrounds from 
the same experiment which attempts to measure the CMB 
anisotropies. 

As already pointed out in the discussion on quality fac- 
tors, polarised foregrounds from dust and synchrotron have 
a very good chance of being extracted as well as the CMB 
signal from future experiments. The basic requirement for 
extracting a process well is that it dominates the total signal 
at least one frequency of the experiment. Both the future ex- 
periments MAP and Planck will have frequency channels 
which fulfill this condition — MAP's two lowest frequency 
channels, centered at 22 and 30 GHz, and LFI's 30 GHz 
channel are likely to be dominated by polarised synchrotron; 
the 545 GHz channel of Planck will act as a template for 
polarised dust. In Figs. ([]) and ( |lo| ) we plot the fractional 
error on the power spectra of various Galactic polarised 
foregrounds. The polarised component of dust can be ex- 
tracted even better than the CMB using HFI because of 
its 545 GHz channel. Similarly, the presence of the lowest 
frequency channel on MAP and LFI should enable one to 
extract the synchrotron signal for I < 100. It should be no- 
ticed that for all the experiments the fractional errors on ET 
cross-correlation are larger than the E power spectrum. It 
is because the ET cross-correlation of foregrounds is mixed 
with a relatively large CMB ET cross-correlation signal and 
therefore does not dominate the total signal as the .E-mode 
foreground polarisation. Another noticeable feature of the 
figures is that MAP performs as well or even better than LFI 
and Planck in extracting the polarised component of syn- 
chrotron. It is because MAP has the lowest frequency chan- 
nel at 22 GHz. Therefore, despite lower sensitivity, MAP has 
better frequency coverage at frequencies which are domi- 
nated by polarised synchrotron signal. 

It should be pointed out that the precision of extraction 
of a given process depends crucially on the other competing 
processes. For instance, the unpolarised component of syn- 
chrotron cannot be extracted as well the unpolarised com- 
ponent not only because the unpolarised component never 
dominates the total signal at any frequency but also because 
it is very difficult to extract it from a comparable level of the 
free-free signal. Therefore, the quality of extraction of syn- 
chrotron polarisation shall depend quite sensitively on the 
presence of other sources of polarised foreground like po- 
larised free-free emission or the extragalactic radio sources 
which are seen to be polarised at 5; 20% level (Saikia & 
Salter 1988). The dust polarised signal, however, is unlikely 
to be affected as there are no known sources of polarised 
foregrounds in the high frequency range covered by HFI, 
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except infra-red point sources which might be polarised but 
will affect the signal only for I > 1000 (Toffolati et al. 1997). 



4.3 Detection of B-mode polarisation 

An unambiguous way to infer the presence of gravitational 
waves in the early universe is through the detection of B- 
mode polarisation (Seljak 1997, Kamionkowski & Kosowsky 
1997, Seljak & Zaldarriaga 1997). (Gravitational lensing can 
generate a B-mode signal from a purely B-mode primary in- 
put . However, this signal is weak and peaks around I ~ 1000 
(Zaldarriaga & Seljak 1998, Bernardeau 1998). We neglect it 
in this paper.) This signal is much smaller than the B-mode 
polarisation and is negligible for £ J; 100, but it is poten- 
tially detectable because of the high sensitivity of Planck. 
We plot the fractional errors on this quantity in Fig. ( ll|) for 
the specifications of Planck. The base model is tilted CDM 
with the scalar index n s = 0.9, the tensor index n t = —0.1, 
and ratio of tensor to scalar quadrapole T/S = 0.7. This 
signal is barely above the noise level and is comparable to 
the foreground contamination (Prunet et al. 1998). How- 
ever, Fig. ([ll]) suggests that it might be possible to get a 
marginal detection with Planck at low £. 

A noteworthy feature of the B-mode detection is that 
the Wiener filtering result is quite comparable to the perfor- 
mance of the combined sensitivity of all the channels. This 
is largely because the B-mode signal does not correlate with 
any other signal and is therefore free of errors coming from 
cross-correlations with foregrounds of other fields (T and B). 
Our results suggest that we expect the answer to lie between 
the case of best channel and the combined sensitivity of all 
the channels. With band-power estimates with a 20% loga- 
rithmic band, the errors drop to < 25% for 20 < £ ~ 100. 
These are shown in Fig. (112). 

Also, as seen in Fig. ( O ) , the B-mode dust polarisation 
can be extracted with much better precision than the CMB 
signal. It is because though the CMB B-mode signal is more 
than an order of magnitude below the B-mode signal, we 
took them to be comparable for foregrounds. This would 
mean that the B-mode foregrounds will dominate the signal 
at most of the frequencies of the future experiments, and 
therefore can be extracted better than the CMB signal. The 
synchrotron B-mode signal, however, is seen to be much 
harder to extract. 



5 CONCLUSION AND DISCUSSION 

We devised a multi-frequency Wiener filtering method to 
consider the effect of Galactic polarised foregrounds on 
the detection of CMB B-mode polarisation and ET cross- 
correlation using future CMB missions. Our results can be 
summarized as: 

(i) The foregrounds can be subtracted well enough for 
the LFI and HFI aboard Planck to detect the E and ET 
signal with signal-to-noise ~ 2-10 for most of the Grange in 
50 < I <, 1000. 

(ii) The foregrounds are likely to greatly undermine the 
performance of MAP. It seems unlikely that MAP could de- 
tect either E or ET signal for individual modes. However, 



by taking band-power estimates with a 20% logarithmic in- 
terval, noise levels reduce sufficiently to allow a marginal 
detection of the ET signal near I ~ 300. 

(iii) The power spectra of B-mode polarised dust can be 
extracted for I ^ 1000 range using Planck HFI with signal- 
to-noise between 1 and 10. The ET cross-correlation of this 
contaminant can be detected with signal-to-noise < 2 for 
100 < I < 1000. The B-mode power spectrum of polarised 
synchrotron, on the other hand, can only be determined by 
either MAP and Planck LFI for t < 100. This suggests 
that both MAP and Planck have a fairly good chance of 
determining the polarised foregrounds which are expected to 
hamper their performances in the detection of the CMB po- 
larisation, at least for a small range of £. 

(iv) The B-mode polarisation, which unambiguously es- 
tablishes the presence of stochastic gravitational waves from 
the inflationary era, can be detected with signal-to-noise ~ 1 
by Planck for I Ss 100. However, band-power estimates 
with a 20% logarithmic band will enable its detection with 
signal-to-noise ~ 2-4 for 20 < i < 100. 

Another possible way to detect very small signals (B 
and ET signal with MAP and B signal with Planck) is to 
image a small fraction of the sky at high Galactic latitudes 
for longer periods. It is expected that Planck will image 
around 1% of the sky at high Galactic latitudes for periods 
5-6 times the all sky average. One advantage of this ap- 
proach is that foreground level at high Galactic latitudes is 
much smaller than the all sky average (e.g. dust emission is 
smaller by nearly a factor of 10 at Galactic latitudes north of 
70° as compared to the all sky average we use in this paper). 
However, this will increase the covariance of power spectra 
by a factor of U / s k y . In light of our results one could ask 
whether it is preferable to deep-image a part of the sky with 
small foregrounds or one should decrease the covariance by 
covering a larger fraction of the sky, albeit with higher fore- 
ground levels. It is seen from Figs g and ^ that the expected 
errors on the power spectra of ET and E are, for £ ~ 100, a 
factor of 6 more than the best performance of MAP without 
foregrounds. It might be possible in this case to gain at least 
a factor of 2 (in addition to sensitivity gained by decrease 
in pixel noise from longer integration time) by imaging 10% 
of the sky at higher Galactic latitudes where the foreground 
level is expected to be much smaller. However, this is not the 
case for the detection of B-mode polarisation by Planck. 
As seen in Fig [H], Wiener filtering extracts the signal almost 
as well as one could get using the combined sensitivity of all 
the channels. This means that if 1% of the sky is integrated 
for a period ~ 5 times more than the all sky average, the 
fractional error will only increase by a factor of ~ 4 because 
a decrease in foreground level will not affect the fractional 
errors. 

In going from Eqs. @ to Eqs. @, @, and we 
assumed the foregrounds to be Gaussian, i.e., the 4-point 
correlation function was assumed to be expressible as a com- 
bination of 2-point functions only. This assumption is erro- 
neous for the foregrounds. The irreducible 4-point function 
from foregrounds can add to the covariances thereby enhanc- 
ing the errors. We cannot quantify this increase within the 
framework of this work. However, as pointed out above some 
of the frequency channels in the future experiments will be 
dominated by foregrounds. A direct analysis of these maps 
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is likely to reveal the non-Gaussian nature of foregrounds, 
which could then be used to quantify errors on CMB extrac- 
tion. 

Dependence on input model : The results presented in 
the paper obviously depend on the choice of model chosen 
for generating the CMB fluctuations. We used a variant of 
CDM model with reionization for the E and ET extractions. 
Reionization suppresses anisotropies at small angular scales 
by exp(— 2r), which is not a major effect for r = 0.1 that 
we choose. The most important effect of reionization is to 
generate new polarisation anisotropies at £ ~ 2-30 depend- 
ing on the optical depth. In the absence of reionization, the 
signal will be smaller by more than an order of magnitude 
for I ^ 15. Therefore, the validity of our results will be very 
sensitive to the underlying model for £ < 15. 

The detection of polarisation anisotropies at smaller 
scales should not depend so much on the choice of model, be- 
cause most variants of sCDM models give comparable level 
of these anisotropies for £ > 20. It should be true unless the 
optical depth to the large scattering surface is large (r ~ 1) 
or that the large scale anisotropies are dominated by ten- 
sor anisotropies. However, the current data for temperature 
anisotropies already suggest that the first Doppler peak is 
even higher than predicted by sCDM model (Netterfield et 
al. 1997), which rules out the possibility of strong reion- 
ization and makes it difficult for the scalar anisotropies to 
be sub-dominant to contribution from gravitational waves. 
Therefore, it is safe to conclude that our predictions for 
£ > 20 will not be seriously affected by a change in the 
underlying model. 

For B-mode polarisation we considered a model with 
tensor to scalar quadruple ratio T/S — 0.7. As mentioned 
above, this will lower the contribution from scalar modes 
and therefore will reduce the signal-to-noise for the detec- 
tion of _E-mode polarisation. However, the B-mode signal is 
roughly proportional to the value of T/S within the frame- 
work of inflationary models which require T/S ~ —7nr- As 
shown above, the B-mode foregrounds can be subtracted 
quite efficiently and for this model the signal-to-noise for 
CMB B-mode signal is ~ 2-4. Therefore, it might be possi- 
ble to detect a signal with a value as small as T/S ~ 0.2-0.3. 

The Wiener filtering method assumes a priori knowl- 
edge of the power spectra of CMB and foregrounds as well 
as the frequency dependence of foregrounds. Therefore, the 
error of extraction using this method does not include the er- 
ror in evaluating these input quantities. Future experiments 
will make multi-frequency maps with millions of pixels. This 
will make it difficult to apply the usual maximum likelihood 
technique to extract various power spectra. Fast methods 
for tackling this problem are currently being developed (Oh 
et al. 1998). In future, it should become feasible to quantify 
errors on the extraction of the power spectra from high res- 
olution multiple frequency maps. It should then be possible 
to revise our estimates of the errors. 
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Figure 1. First Panel : The power spectrum of the expected E- 
mode CMB signal is shown along with the dust and synchrotron 
polarised foregrounds at different frequencies. The dotted line 
show the CMB B-mode power spectra. The dashed lines corre- 
spond to dust polarised emission (Eq. ( pq ) while the dot-dashed 
curves correspond to polarised synchrotron (Eq. (^)). The four 
lines for foregrounds give the expected level of contamination 
at 44 GHz, 100 GHz, 143 GHz and 217 GHz. The corresponding 
curves at these frequencies go from top to bottom (bottom to 
top) for synchrotron (dust) as the frequency is increased. The 
thick solid line gives the temperature power spectrum for the 
same underlying model, which is taken to be sCDM with r = 0.1. 
Second Panel : Same as the First Panel for ET cross-correlations. 



Planck - HFI 



Planck - l.FI + HFI 



Figure 2. Quality factors for the temperature fluctuations. Solid, 
dashed, and dot-dashed lines stand for CMB, dust and syn- 
chrotron respectively. As can be seen by comparing the effective 
windows in the case of different experiments, despite its lower 
sensitivity MAP does as well as the LFI at recovering the syn- 
chrotron emission because of its larger frequency coverage; but 
since this emission is relatively weak, the LFI CMB transmission 
is better. Also note the improvement on the foreground recovery 
when the LFI and HFI are combined, although the CMB recovery 
is not improved by adding the LFI data. 
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Figure 3. Quality factors for the _E-mode polarisation. Solid, 
dashed, and dot-dashed lines stand for CMB, dust and syn- 
chrotron respectively. Note that all the experiments extract the 
foregrounds better than the CMB for small I. Planck HFI ex- 
tracts the polarised dust better than the CMB at almost all I. 
This suggest that polarised foregrounds arc likely to be extracted 
as well or better than the CMB. Though this is partly due of our 
assumed level of foregrounds, this feature reflects the frequency 
coverage of these experiments 



Figure 4. Relative fractional error on temperature power spec- 
tra. The first panel shows the expected performance of the full 
Planck mission in extracting this signal. The Solid line illus- 
trates the Wiener filtering case, i. e. the foregrounds are included 
and subtracted using the Wiener filtering method described in 
the text. The dotted line corresponds to the case if the signal 
is extracted using only the best channel of Planck (143 GHz) 
neglecting all the foregrounds. The other panels show the perfor- 
mance of other experiments relative to the Planck- Wiener case 
(the solid line in the first panel). In addition to the Wiener and 
best channel case (the same line style as the first panel), we also 
show the errors using the combined sensitivity of all the channels 
{dashed line) for various e: 
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Figure 5. Fractional errors in extracting the E mode polarisa- 
tion. The solid, dotted, and dashed lines correspond to the ex- 
pected performances of Wiener filtering, best channel, and com- 
bined sensitivity of all the channels of a given experiment, respec- 
tively. 
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Figure 6. Same as Fig. g for the ET cross-correlation. The sharp 
spikes in the figure merely indicate the values of £ at which the 
signal vanishes. 



© 0000 RAS, MNRAS 000, 000-000 



12 Bouchet et al. 



Dust extrr.etrcn 




Figure 7. £?-modc band power estimates. We plot the _E-mode 
signal for the underlying model {solid line) and show the expected 
la measurements [dashed , dotted and dot-dashed line for HFI, 
LFI and MAP respectively) when band power with a logarithmic 
interval A£/£ = 0.2 is taken. 





Figure 9. Fractional errors on the extraction of dust E and ET 
power spectra with HFI and full Planck . The lower two curves 
correspond to E power spectrum {solid and dotted lines are for 
HFI and Planck respectively) while the upper curves give the 
errors for ET cross-correlation, with the dashed line for Planck 
and the dot-dashed line for Planck HFI. 

Synchrotron extraction 




Synchrotron extraction 




Figure 10. Extraction of synchrotron E (Left Panel) and ET 
(Right Panel) power spectra with LFI {dotted line), full Planck 
{solid line), and MAP {dashed line). 



Figure 8. Same as Fig. [?| for the ET cross-correlation. 
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Figure 11. Fractional errors for B-mode polarisation arc shown 
in the first three panels. The solid , dotted , and dashed curves 
correspond to Wiener, best channel, and combined sensitivity, 
respectively. The forth panel shows the corresponding errors in 
the extraction of the S-mode component of the polarised dust 
(solid line) and synchrotron (dotted line). The underlying model 
was taken to be a CDM model with n s = 0.9, T/S = 0.7, and 
tit = 0.1 
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